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A systematic approach is presented to construct non-homogeneous two- and three- 
qubit Bell-type inequalities. When projector-like terms are subtracted from homoge- 
neous two-qubit CHSH polynomial, non-homogeneous inequalities are attained and 
Qh| the maximal quantum mechanical violation asymptotically equals a constant with 
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the subtracted terms becoming sufficiently large. In the case of three-qubit system, it 
is found that most significant three-qubit inequalities presented in literature can be 
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recovered in our framework. We aslo discuss the behavior of such inequalities in the 



• loophole-free Bell test and obtain corresponding thresholds of detection efficiency. 
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I. INTRODUCTION 



Bell's original inequality reveals the conflict of quantum mechanics and local hidden 
variable (LHV) theory [1] . Since then, various forms of Bell-type inequalities have been de- 
rived. The most well-known is Clauser-Horne-Shimony-Holt (CHSH) inequality [2]. Mermin- 
Ardehali-Belinskii-Klyshko (MABK) generalized the result to the case of A" qubits [3]. Sub- 
sequently, Werner and Zukowski (WWZB) presented explicit description of A^-qubit Bell- 
correlation inequalities for two dichotomic observables per site [4, 5]. It should be stressed 
that all these results only concern with the "full" correlation function, that is, the expecta- 
tion value of the product of all A" local observables. We call this class of Bell-type inequalities 
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as homogeneous inequalities. 

Besides full correlations, one can take into account partial correlations, which involve 
the product of not N but n < N observables. It is meaningful to construct such non- 
homogeneous Bell-type inequalities that incorporate partial correlations as well as full cor- 
relations. The main reason is as follows. Gisin proved that the bipartite pure entangled 
states violate CHSH inequality [6, 7]. The problem whether Gisin theorem can be general- 
ized for an arbitrary n-partite pure entangled states remains open. In the three-qubit case 
there are generalized GHZ states that do not violate the MABK inequalities [8, 9]. More 
generally it has been shown that these states do not violate any Bell-type inequality for 
n-partite correlation functions (that is, homogeneous inequality) for experiments involving 
two dichotomic observables per site [9]. Subsequently some Bell-type inequalities, which 
bear the non-homogeneous form, were presented and can be violated by generalized GHZ 
states [10, 11]. Besides, non-homogeneous inequalities have appeared in the discussion of 
the non-locality of cluster states [12, 13]. However one can hardly make clear how and 
why these inequalities come out. In other words, there is no systematic formulation of such 
inequalities. 

Motivated by this issue, we present in this paper a feasible method to construct non- 
homogeneous Bell-type inequalities. The key point is so simple that we can describe it very 
briefly in one sentence, that is, non- homogeneous inequalities are derived by subtracting 
some projector-like terms from the homogeneous ones. 

We will first present the construction of two-qubit non-homogeneous inequalities. The 
newly obtained inequalities are weaker than CHSH inequality in that the quantum mechan- 
ical violation of the former is smaller than that of the latter. Nevertheless, they pave the 
way for the considering more general cases. Consequently we attain so many three-qubit 
non- homogeneous inequalities and some meaningful results. In Ref. [11], Pitowsky and 
Svozil have presented several optimal three-qubit Bell-type inequalities in the sense that 
they represent the best possible upper bounds for the conceivable classical probabilities. 
These inequalities are filtrated from tens of thousands of inequalities describing the facets of 
classical correlation polytope. Though these inequalities appear very intricate, they can be 
constructed systematically in our framework. Additionally the two Bell-type inequalities, 
which were presented in Ref. [11] to disclose the non-locality of generalized three-qubit GHZ 
states, are included in our results. 
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When these inequalities are utilized to display the non-locality experimentally, one has to 
consider the effect of non-ideal detector. Using the software Mathematica and MatLab, we 
numerically analyze the detection efficiency limit of the above-mentioned three-qubit Bell- 
type inequalities. Three detectors may have the same efficiency or not [15]. With respect to 
these two different cases, we attain the thresholds of detection efficiency Some inequalities 
have the advantage that the efficiency of one detector can be arbitrarily low if that of the 
other two detectors satisfy certain conditions. 

II. TWO-QUBIT NON-HOMOGENEOUS BELL-TYPE INEQUALITIES 

We start by briefly recalling LHV model and Bell-type inequalities. Bell-type inequalities 
always refer to correlations between two or more sites. In the two-qubit case, each of 
two space-separated observers, A and B, gets a qubit and measures two ±l-valued local 
observables, denoted A\, A 2 for observer A, and B\, B 2 for observer B. The outcomes of 
measurement are labeled by = ±1 and bj = ±1 for i, j = 1, 2. 

A typical joint probability can be expressed as P(a 2 , bi\A 2 , Bi), where after the vertical 
bar we write the observables chosen at two sites, and before the bar the particular outcomes. 

In the formalism of local hidden variable (LHV) theory, there exists a hidden variable A 
which takes values in space A. With the presence of A, the probability of measuring A4 and 
obtaining the outcome a, is represented by P(ai\Ai, A). Similarly for P(bj\Bj, A). Given A, 
one can calculate the mean value of aj(A) from the probability P(ai\Ai, A), 

a l (A)=P(+l|A i ,A)-P(-l|A i ,A). (1) 

Similarly for bj(A). 

Furthermore LHV theory requires that with the presence of A the joint probability 
P(aj, bj\Ai, Bj, A) is factorisable, that is, 

P{a i ,b j \A i ,B j ,X)=P{a i \A i ,X)P{b j \B j ,X), (2) 

where P(ai\Ai, A) is independent of observable Bj and its outcome bj, and P(bj\Bj, A) inde- 
pendent of Ai and a;. 

Given the probability measure /i on A, using (2), one can compute the joint probability 
P(a i ,b j \A i ,B j ) = [ P(a i \A i ,X)P(b j \B j ,X)n(X)d\. (3) 

J A 
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Consequently the expectation value of joint measurement is 

3i,bj 

= [ J2a i b j P(a i \A i ,\)P(b j \B j ,\)ii(\)d\ 



a i (A)b j (A)/x(A)dA. (4) 

A 

In the context of hidden variable model, correlation function is defined as 

E(X) = ai(A)bi(A) - a 1 (A)b 2 (A) 

-a 2 (A)b 1 (A)-a 2 (A)b 2 (A). (5) 

One can see that |-E(A)| < 2, and the inequality is saturated only if a"j(A) and bj(A) take the 
extremal values of ±1. CHSH inequality is attained by averaging over hidden variable A. 

| (A 1 B 1 ) lhv - (A 1 B 2 ) lhv - (A 2 B 1 ) lhy - (A 2 B 2 ) lhv \ < 2. (6) 

In the following discussion, it is convenient to adopt |-E(A)| < 2 to represent CHSH 
inequality, that is, 

- 2 < a±bi - aib 2 - a 2 bi - a 2 b 2 < 2, (7) 

where we replace a"j(A) and bj(X) with ai and bj for simplicity. We also make the following 
definitions. 

E l = ai&i - a x b 2 - a 2 b x - a 2 b 2 , (8) 

E 2 = -ai&i + a x b 2 - a 2 b x - a 2 b 2 , (9) 

E 3 = -cq&i - a x b 2 + a 2 b x - a 2 b 2 , (10) 

E 4 = -ai&i - a\b 2 - a 2 b x + a 2 b 2 . (11) 

Each E k is called (homogeneous) CHSH polynomial and has the same bound of 2, i.e., 
\E k \ < 2 for k = 1,2,3,4. 

A. Method 



Now we start to construct non-homogeneous polynomial, denoted E '. We require that 
E' satisfies two conditions: 
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(i) E' < 2. That means in the LHV upper bound of E' is 2. 

(ii) In the case of quantum mechanics, E' is replaced with the operator form £', and the 
expectation value given by (^') qm can be larger than 2. 

In other words, E' < 2 is a Bell-type inequality. 

The key point of the method presented here is so transparent that at first sight it seems 
to be trivial: since every CHSH polynomial satisfies < 2, it is evident that E^ minus 
some non-negative terms must be equal to or less than 2. These non-negative terms are 
chosen to be 

Pt = ~, Ft = -• 12 

(ii 2 2 

A typical non-homogeneous polynomial is then given by 

E'(r)=E x -rP+P+, (13) 

where r is any non-negative number. 

Considering condition (i), we note that E' < 2 and E' reaches the maximal value of 2 at 
specific extremal values of <2j and bj. For example, when a x — 1, 02 = 1, b\ — 1 and b 2 = — 1, 
E' equals 2. So (13) satisfies condition (i). 

We now show that condition (ii) can also be fulfilled. In quantum-mechanical case, and 
bj in (12) are replaced with spin observables, namely, =^ Ai, bj =3- Bj. Correspondingly, 
we have 

E 1 =>£ x = A X B X - A X B 2 - A 2 B t - A 2 B 2 , (15) 

where / is 2 x 2 identity matrix. 

Then the quantum mechanical form of (13) is 

£'{r)=£ x -rVZ 2 V+, (16) 

Any two-qubit entangled pure state can be expressed as the form of Schmidt decomposition 
(see, for example, [16]), 

= cos£|00) + sinf|ll), fe(0,7r/2). (17) 
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We let 

A x = a x , A 2 = —a x , 
Bi = a x sin 9 + a z cos 9, 
B2 = —o x sin 9 + a z cos 9. 
The expectation value of £' is given by 



(18) 



= (2 + ^sin 2 ^) costf 



+ 2 sin 2£ sin 9 — — sin £. 

When = arctan ^rsiSg ' Wl^> 

acquires the maximal value 



(19) 



r . 2 x2 



2 + - sin z £ ) +4 sin 2 2£ 



1/2 



2^' 



(20) 



which is obviously larger than 2. Therefore (13) represents a series of non-homogeneous 
inequalities and can be used to detect non-locality 

The above procedure can be formulated in more general forms. Let's consider the follow- 
ing non-homogeneous polynomial. 



E*(8,t) = E A -sP+P+-tP+P+ t 



(21) 



where s and t are non-negative numbers. This time we select E4 to construct E". It is just 
for the later convenience and not necessary. Obviously we have E" < 2 and the bound is 
attained for some extremal points, say, a\ — a-i — 1 and 61 = b 2 = — 1- In order to show 
the violation of E" < 2 in quantum mechanics, we evaluate the expectation value 
where |\&) is given by (17) and £" is the quantum conterpart of E" . 
For £ e [0, vr/4], we choose the observables as 



Ax = a z , A 2 



Bi = o x sin 9 + o z cos 9, 
B 2 = —<J X sin 9 + a z cos 9. 



(22) 



For this choice, we have 



+ 



t + 



[s + i) cos2£] cos9 
1 



j sin 2£ sin 6* — — (s + i + scos2£) 



(23) 
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The maximum of over all 9 is given by 

/ 1 (s,t,0 = max(*|£"|*) 

all 

= -^(s + t + scos 2f) 

+ J [ (8 + s + (s + t) cos 20 2 + (8 + t) 2 sin 2 2f ] 1/2 ■ (24) 
It can be found that if s and t satisfy 

32 - st + (32 - st + 8t) cos 2£ > 0, £ e [0, tt/4] , (25) 

the function /i(s, i, £) will be larger than 2 (when £ = 0, fa = 2). 
For £ e [7r/4, 7r/2] , we choose 



A\ = —a z , A 2 = —a x , 
Bi = <j x sin 6 + a z cos 9, 
B2 = —<J X sin 6 + a z cos 

Similarly, we have 

fa(s,t,Z)=max(*\S"\V) 

all y 

= — - (s + 1 — s cos 2£ ) 



If s and t satisfy 



(26) 



+ J [(8 + s - (s + t) cos 2£) 2 + (8 + t) 2 sin 2 2f ] 1/2 . (27) 



32 - st - (32 - st + 8t) cos 2£ > 0, £ G [tt/4, tt/2] , (28) 



then fa(s,t,£) is larger than 2 (when = 7r/2, / 2 = 2). 

It can be seen from (25) and (28) that both of /1 and fa are larger than 2 if st < 32. So, 
under this condition, {£") violate LHV upper bound. 

Some remarks are needed here. In the above calculation, we do not intend to find the 
maximal violation of LHV bound. The choices of observable, given by (18), (22) and (26), do 
not cover all possible local measurements. They are so selected as to simplify the calculation. 
Therefore (\I / |^ / |\I / ) max and ( i f?\£"\ty) raBX do not represent the maximal violation of LHV 
bound. In the same way, the condition of st < 32 is just required in the case considered. 
In the following subsection, we discuss the maximal violation of inequalities E' < 2 and 
E" < 2. 
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B. Violations of inequality 

Two typical forms of non-homogeneous inequality have been presented in (13) and (21). 
We now analyze the violation of them in quantum mechanics. The quantum bound can be 
computed by means of Min-Max principle [17], which states that for self-adjoin transfor- 
mations the operator norm is bounded by the minimal and maximal eigenvalues. For the 
purpose of clarity, we choose the directions of local measurement to be orthogonal. In other 
words, A\ and A 2 anti-commute. So do B\ and B 2 . Furthermore, we let 

A 1 = B 1 = a x , A 2 = B 2 = a y . (29) 

The operator S'(r), which is quantum correspondent of E'(r), is written as 

/-l - r -1 -1 1 \ 

S'(r) = 



(30) 



-1 111 
-1111 

\ 1 11-1/ 

The eigenvalue of £'(r), denoted s(r), is determined by the equation 

e 4 + re 3 - (r + 8)e 2 - Are = 0. (31) 

Among the solutions of (31), the largest one must be between 2 and 3. To see this, note that 
the left-hand side of (31) is negative when e = 2, while positive and monotone increasing 
when e > 3. 

The numerical result of the largest eigenvalue of S'(r) is plotted in Fig. 1. It can be 
seen that the largest eigenvalue asymptotically equals a constant (about 2.5) for sufficiently 
large r. The asymptotic value can be calculated from (31). In fact, when r — > 00, the terms 
containing r should be finite. It follows that 

e 3 - e 2 - Ae = 0. (32) 

The solution ~ 2.56 is what we want. 

Similar analysis shows that the largest eigenvalue of S"(s,t) is also between 2 and 3, and 
asymptotically equals 2.43 when s,t —>■ 00. 

In a word, the inequalities E'(r) < 2 and E"(s,t) < 2 can be violated by quantum 
mechanics for arbitrary non-negative numbers of (r, s, t). The maximal violation corresponds 
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FIG. 1: (Color online) The maximal eigenvalue of S'{r). When r tends to infinity, the asymptotic 
value is ±±^S ~ 2.56. 

to the case of r = s = t = 0. With r, s and t sufficiently large, the violation is close to a 
constant. 

III. THREE-QUBIT NON-HOMOGENEOUS INEQUALITIES 

Based on the results obtained previously, we will in this section construct three-qubit 
non- homogeneous inequalities. 

A. General results 

We now introduce the third qubit C. Two observables C\ and C<i are also ±I-valued. Just 
like Eq. (1), we define 



The general form of three-qubit non-homogeneous polynomial can be expressed as 



c i (A) = P(+l|C i ,A)-P(-l|C i ,A). 



(33) 



F(a, 6)ci + G(a, b)c 2 + H(a, b), 



(34) 
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where F, G and H are polynomials associated with qubit A and B. We will construct three- 
qubit inequalities by the following steps. Firstly, we set the upper bound of (34) to be an 
arbitrary number, for example, 

F(a, 6)ci + G(a, b)c 2 + H(a, b) < 2 + u, (35) 

where we let u > simply for the technical simplicity. Secondly we let C\ and c 2 take 
extremal values of ±1. (35) will reduce to four inequalities with the upper bound of 2 + u. 
We let three of them be the known two-qubit Bell-type inequalities (homogeneous or non- 
homogeneous). For example, we make the following choice. 

F + G + H = E" + u<2 + u, (36) 

F-G + H = E' + u<2 + u } (37) 

-F-G + H = ^p(-£ 4 ) < 2 + u. (38) 

So F, G and H can be expressed in terms of E'(r), E"(s,t) and £4. 

„ 2 + ta E' u 

F -— E ' + Y + r (39) 

G = (40) 

ff = __£,+ _ + _. (41) 

Finally, there remains a reduced inequality, namely, — F + G + H<2 + u. We use it to 
determine the possible values of r, s and t. From (39), (40) and (41), it follows that 

2-1-7/ 

-F + G + H = -—-^E± + E"-E'. (42) 

The requirement that the upper bound of (42) is 2 + u leads to the following restrictions on 
r, s and t. 

r < 4 + 2u, (43a) 

r - s < 2u, (43b) 

r-t<2u, (43c) 

r - s - t < 0. (43d) 

In a word, we require that four reduced two-qubit polynomials, namely ±F ± G + H, have 
the same upper bound of 2 + u. We point out that the choice given by (36)-(38) is not 
unique. What is more, some reduced inequality are not necessarily required to be Bell-type 
inequalities, as will be seen in the next subsection. 
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B. Specific inequalities 

We present some examples to illustrate the application of our method. 

1. u = 0, r = s = t = 

It follows that 

E' — Ex, E" = E 4 , (44) 
F = ~(£ 4 + Ex), G = 1 -{E, - Ex), H = 0. (45) 

Then we have 

- aibic 2 - GH&2C1 - a 2 bic 1 + a 2 b 2 c 2 < 2. (46) 

It is the well-known three-qubit MABK inequality. The largest eigenvalue of the quantum 
correspondent is 4, which means the violation factor is | = 2. 

2. u = 2,r = 8,s = t = A 

In this case, conditions (43) are satisfied and saturated. Direct calculation gives 

— at — a 2 — bi — b 2 + ai&i — 02^2 — a\c 2 

- 2a 2 d + a 2 c 2 - bic 2 - 2b 2 C\ + &2C2 

— ai&iCi — 2a\b\c 2 — "ia\b 2 c\ — a\b 2 c 2 

- 3a 2 biCi - a 2 bic 2 - a 2 b 2 C\ + Aa 2 b 2 c 2 < 8. (47) 
In Ref. [14], Pitowsky and Svozil enumerate some new Bell-type inequalities. In their 

notation, one of the inequalities (Eq. (5) in [14]) is 

- P(Ax) - 2P(B 1 ) - 2P(C 1 ) + 2P(A h B x ) 

+ 2P(A U d) + P(Ax, B 2 ) + P(A U C 2 ) + P(A 2 , Bx) 
+ P(A 2 , Cx) - P(A 2 , B 2 ) - P(A 2 , C 2 ) + 2P(Bx, Cx) 
+ 2P(B 2 , Cx) + 2P(Bx, C 2 ) - 2P(B 2 , C 2 ) 

- P(Ax, Bx, Cx) - 2P{A 2 , Bx, Cx) - 3P(A U B 2 , Cx) 

- 3P(Ax, Bx, C 2 ) - P(A 2 , B 2 , Cx) - P(A 2 , Bx, C 2 ) 
-P(Ax,B 2 ,C 2 )+AP(A 2 ,B 2 ,C 2 ) <0. (48) 
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Here P(Ax) denotes the probability that observer A measures observable A\ and obtains the 
outcome +1, that is, in our notation, P(ai = +l|Ai). We know that 

P(ai = +l|Ai)= / P(ai = +l|Ai,A)fz(A)dA 

i±^.^(A)dA. (49) 
Similarly P(Ai, P>i) refers to P(a\ = +1, bi = +1 1 Ai, Bi). And we have 

P(ai = +l,bi = +l|Ai,Bi) 



A 



: / P(a 1 = +l,b 1 =+l\A 1 ,B 1 ,X)ii{\)d\ 

J A 

[ P(a 1 = +l\A 1 ,\)P(b 1 = +l\B 1 ,\)fi(\)dX 

J A 



A 1+ 2 "' 1+ 2 6 V (A)rfA. (50) 

Then the inequality (48) can be rewritten in terms of a^, 6j, c^. 

We find that, under permutation b — > a, a — > c and c — > b, Pitowsky's inequality (48) 
can be transformed to (47). Furthermore, the other inequalities presented in [ 1] can be 
recovered in our framework. 

In quantum mechanics, the largest eigenvalue of the operator corresponding to the left- 
hand side of (47) is about 12.87 (numerically), which means the violation factor is about 
w 1.61. 

We also evaluate various violation factors for different values of (u,r,s,t). Numerical 
results reveal such phenomena: (i) When u < 2, the optimal inequality, which can give rise 
to the maximal violation of LHV bound, corresponds to the case of r = 4u and s = t — 2u, 
that is, conditions (43b), (43c) and (43d) are saturated, (ii) when u > 2, the optimal 
inequality is attained for r = 4 + 2u and s = t = 2 + u, that is, conditions (43a) and (43d) 
are saturated. In Fig. 2, we plot the maximal violation factor for different values of u. It 
can be seen that the violation factor tends to be a constant (about 1.27) for large u. 
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FIG. 2: (Color online) The maximal violation factor for different u. When u sufficiently large, the 
asymptotic value is about 1.27. 



3. Another case 



Let's recall the steps by which the three-qubit non-homogeneous inequalities are estab- 
lished. The key point is to let three reduced polynomials be related to three two-qubit 
Bell- type inequalities, which is demonstrated in Eqns. (36), (37) and (38). This require- 
ment can be relaxed. For example, we consider the following procedure. 

F + G + H = 2 - P+P+ - P+F+ < 2, (51) 
-F + G + H = E 4 <2, (52) 
-F-G + H = 2-P-P b --P-p-<2. (53) 

The LHV upper bounds of above polynomials are the same. But only (52) is a real Bell-type 
inequality. From these equations, we get the following inequality. 

— di&i — d\b 2 — a 2 bi — a 2 b 2 — a\C\ 

— aic 2 — a 2 c\ — a 2 c 2 — b\C\ — b\c 2 

- b 2 ci - b 2 c 2 + aibiCi - a\b 2 c 2 

- a 2 bic 2 - a 2 b 2 c 2 + 2a 2 b 2 c 2 <4 , KA s 
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This inequality was originally presented in Ref . [1 1] and said to be violated by all three-qubit 
pure entangled states (see Eq. (6) in [11]). 

IV. THRESHOLDS OF DETECTION EFFICIENCY 

In the actual experiment performed on Bell-type test, there are mainly two kinds of 
loophole to overcome, that is, locality loophole and detection loophole [18]. Locality loophole 
has been closed in several Bell experiments [19]. However, detection loophole is difficult to 
overcome so far [: 0]. In order to perform a loophole- free Bell test, one must make clear 
the threshold detection efficiency. Many useful results have been obtained for different Bell 
inequalities in two-qubit system [15, 21]. In this section we discuss the thresholds of detect 
efficiency for six three-qubit inequalities. Four of them come from Ref. [J I] and are labeled 
from PI-(2) to PI-(5) where the number denotes the equation number in the reference. The 
other two inequalities, labeled CI-(2) and CI- (6), are respectively Eq. (2) and Eq. (6) in 
Ref. [11]. 

We take inequality (48), i.e., PI- (5), as an illustrative example. As said before, P(Ai) 
is the probability of finding the outcome ai = +1, and P(A\,Bi) is a joint probability of 
finding the outcomes ai = +1 and t>i = +1. Other probabilities have the similar meaning 
and all of them are defined in ideal situation. 

Now we assume that the detection efficiencies at site A, B and C are rji, r] 2 and 773 re- 
spectively. Then P{A\) should be multiplied by a factor 771. The joint probability P(A\, E>i) 
is replaced by r]ir] 2 P(A 1 , Bi), and P(Ai, B 1: Ci) by 771V72I3 P(Ai, B x , Ci). Similarly for other 
probabilities. Thus we get the modified form of PI- (5), which is utilized to deal with non- 
ideal cases and should be satisfied by LHV model. 

In quantum-mechanical formalism, the probability P(Aj) is replaced by (^f\V^ ® I <S> 
J|#) = (*|^ for some three-qubit state |\&), and the joint probability P{A%, B{) 
by (ty\V+ ®V^®I\^) = (^l^^ <g> ^^y 1 ® I\#). Similarly for other probabilities. Thus we 
obtain the quantum mechanical expression of the modified PI- (5), that is, 

(*|c7|*> < (55) 

where J contains such operators as rjiP^, VityVaVj^, etc. Once we find ($1^71^) > 0, a 
loophole-free experiment can be performed. In other words, we hope to find that the maximal 
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eigenvalue of the operator J is positive. Before we proceed to numerical calculation, it is 
necessary to lower the number of the parameters in J . Since we are to find the eigenvalues 
of J ', local observables can be set to be in xy-plane. Without loss of generality, we let 

Ax — Bx — Cx — a x , 
A2 = cr x cos 9a + o y sin 9 a, 
B2 = <? x cos 9b + o y sin6>£, 
C2 = o'x cos 9 C + cr y sin 9 C . 

Thus the 8x8 matrix J contains six parameters: 771, 772, 773, 9a, 9b, 9c- Then under the 
condition that the maximal eigenvalue of J is positive, we find the threshold of detection 
efficiency numerically by means of Mathematica and MatLab. In the following we discuss 
two cases and give corresponding results. 

a. Symmetric system - In this case, all three detection efficiencies are the same, i.e., 
f]\ = f}i = V3 = V- we obtain the threshold efficiency rf h = 66.8% on which detection 
loophole can be closed. 

b. Asymmetric system - In a asymmetric system, three articles may be detected with 
different probabilities. We obtain the following results. 

(i) . When one of three detectors is perfect, say, r\\ — 1 and the other two are imperfect 

but have the same efficiency, that is, r\\ — 1, r? 2 = i] 3 = rj, the threshold efficiency is 
found to be r/ th = 50%. 

(ii) . Two detectors are perfect, e.g., rji — 772 = 1, and the third is not. In this case, no 

matter how the efficiency of the third detector is low, the inequality can be always 
violated by appropriate choice of detection orientation. 

(iii) . Assuming that 771 = 1, just as (i), we consider such a problem: With 773 arriving at 

its lower bound, how about the value taken by 772? Calculation results show that the 
lower bound of 773 can be zero while 772 must be larger than 50%. Obviously (ii) is the 
specific case of this result. 

Inequalities PI-(2) and PI- (4) have similar behavior, whereas inequalities PI-(3), CI-(2) 
and CI-(6) do not. For example, for CI- (2) the lower bound of 773 is 81.9% and can not be 
zero, and moreover when 773 = 81.9% the other two detector must be perfect. All results are 
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listed in Table I. Obviously inequality PI-(5) is the most optimal one among them in the 
sense that this inequality can endure very low detection efficiency 



TABLE I: Threshold Detection Efficiency 





symmetric 


asymmetric 




V 


Vi 


V2 


V3 






1 


50% 


50% 


PI-(5) 


66.8% 


1 


1 









1 


60% 


60% 


PI-(4) 


71.5% 


1 


1 









1 


77.2% 


77.2% 


PI-(3) 


87.6% 


1 


1 


34% 






1 


66.8% 


66.8% 


PI-(2) 


87% 


1 


1 









1 


87.4% 


87.4% 


CI-(2) 


91.5% 


1 


1 


76.4% 






1 


90.5% 


90.5% 


CI-(6) 


93.6% 


1 


1 


81.9% 



V. CONCLUSION 



In conclusion, we propose a systematic approach of constructing non-homogeneous Bell- 
type inequalities for two- and three-qubit system. In the two-qubit case, non-homogeneous 
inequality is attained by subtracting positive projector-like terms from CHSH polynomial. 
We find that when the subtracted terms are sufficiently "large" the maximal quantum me- 
chanical violation asymptotically tends to be a constant. 

Three-qubit non-homogeneous inequalities are attained by direct generalization of two- 
qubit ones. Most of significant three-qubit inequalities presented in literature are recovered 
in our framework. The method presented in this paper can be generalized to construct Bell- 
type inequalities for multipartite system. The benefit of non-homogeneous inequalities lies 
in considering not only full correlations but also partial ones, which are needed in discussing 
the non-locality of multipartite system. We conjecture that our method may be used to 
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categorize Bell-inequalities with various and complicated forms. 

Additionally, we analyze numerically the detection efficiency thresholds of the previously 
mentioned three-qubit Bell-type inequalities when they are employed to display non-locality 
in quantum states. Under different situations, that is, three detectors may have the same 
or distinct efficiency, we obtain the thresholds of detection efficiency respectively. For some 
inequalities, we find that the efficiency of one detector can be arbitrarily low as long as 
the other detectors satisfy certain conditions. Numerical results will help us find the most 
optimal one which can endure very low detection efficiency. 
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